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I.  INTRODUCTION 


The  flight  stability  of  liquid-filled,  spin-stabilized  projectiles  has  been  considered  for 
/  a  wide  variety  of  conditions,  Originally,  theories  and  experiments  were  centered  about  the 
case  of  large  Reynolds  number  (Re  =  a2<p/ v,  see  List  of  Symbols)  for  cylindrical  payload 
containers.1’2  Non-steady  effects  must  be  considered  for  practical  applications  and  large 
Reynolds  number,  for  example,  spin-up.3’4’5  However,  if  the  Reynolds  number  is  small 
(Re  <  500),  the  effects  of  unsteady  processes  may  be  neglected.  It  is  then  possible  to 
employ  a  steady  finite- difference  solution  to  the  incompressible  Navier- Stokes  equations 
for  a  fixed  precession  angle  and  steady  rates  of  spin  and  precession.  Pressures  have  been 
measured  under  these  conditions  and  can  be  used  to  validate  numerical  simulations.6  7 

This  report  describes  the  steady-state  finite-difference  method  of  Strikwerda8  that  has 
been  extensively  investigated  by  Nusca9  for  fluid-filled  cylindrical  containers,  and  the 
extension  of  this  code  to  non-cylindrical  axisymmetric  geometries.  In  particular,  the  case 
of  a  cylinder  with  rounded  endcaps  (see  Figure  1)  is  investigated. 

II.  BACKGROUND 


In  this  section  a  steady-state  computational  fluid  dynamics  (CFD)  method  and  an 
analytical  simulation  for  cylindrical  containers,  are  reviewed.  Analytical  and  expeiimental 
work  for  non-cylindrical  containers  is  also  referenced.  In  addition,  the  concepts  of  liquid- 
induced  roll  and  yaw  (side)  moment  coefficients,  by  Murphy,1  are  described. 


Murpkp,  CM.,  “Angular  Motion  of  a  Spinning  Projectile  with  a  Viocon.  Liquid  Payload,"  A RBRL-MR-031 94,  U.S. 
A  rmf  Bath, tic  Retearck  Laboratory,  Aberdeen  Proviny  Cronnd,  Maryland,  Augutt  1  98t.  (AD  A1I8676)  Al  to  AIAA 
Journal  of  Gaidanee,  Control,  and  Dynamict,  Vo I.  6,  .Vo.  4,  pp. 380-386,  July-Auyutt  1983. 

Gerber,  N.  and  Seine y,  R.  “Afomtnl  on  •  Liquid- Filled  Sp mniny  and  Nutating  Projectile:  Solid  Body  Rotation,"  AR8RL- 
TR-OH70,  U.S.  Ballistic  Reiearek  Laboratory,  Aberdeen  Prouing  Ground,  Maryland.  February  1983  (AD  A135333). 

r  3  Murphy,  C.  H.,  “ Moment  Induced  by  a  Liquid  Payload  Dun ny  Spin-Up  W.lAonf  a  Critical  Layer,"  ARBFL-TR  03581, 
L’.S.  Army  Ballietic  Reiearek  Laboratory,  Aberdeen  Proviny  Ground,  Maryland,  Auyuet  1984.  (AD  A14S716)  Alto  AIAA 
Journal  of  Guidance,  Control,  and  Dynamic!,  Vol.  8,  No.  3,  pp. 354-359,  May-Junc  1985. 

Gerber  ;V. ,  "Liquid  Moment  on  a  Filled  Coning  Cylinder  Duung  Spin-Up:  Ad  Hoc  Model,"  ARBRL-  TR-3638,  U.S.  Army 
Ballietic  Retearck  Laboratory,  Aberdeen  Proving  Ground,  Maryland,  December  1984 ■  (AD  150180 ) 

3  D'Amico  W.  P.,  “Flight  Data  on  Liquid-Filled  Shell  for  Spin-Up  Inetabilitiei,"  ARBRL-MR  03334,  U.S.  Army  Ballietic 
Reiearek  Laboratory,  Aberdeen  Proving  Ground,  Maryland,  February  1984.  (AD  139136)  Alio  AIAA  Paper  63-3143,  Augutt 
1  983. 

* Nutca  M.  J.,  D’Amico  W.  P.,  and  Beimi,  W.  G.,  “Preiiurt  Menurementi  in  a  Rapidly  Rotating  and  Coning.  Highly 
Imicom  Fluid,”  ARBRL-MR-0SS3S,  U.S.  Army  Bathetic  Reiearek  Laboratory,  Aberdeen  Proving  Ground,  Maryland  November 
1983.  (AD  A  1 36834) 

Hepner,  D.J.,  Kendall,  T.M.,  Davn,  B.S.,  and  Tenly,  W.Y.,  “ Prenurt  Menurementi  in  a  Liquid-Filled  Cylinder  Uimy  A 
Three- Degree- Of- Freedom  Flight  Simulator ,"  ARBRL-MR-3560,  U.S.  Army  Ballutic  Retearck  Laboratory,  Aberdeen  Proving 
Ground.  Maryland.  December  1986  (AD  A177673).  Alt o,  AIAA  Paper  No.  86-t0t6,  1986. 

Strikwerda,  J.  w.,  and  Nagel,  Y.  M.,  “A  Numerical  Method  for  the  Incompmiible  Navier-Stoket  Equationi  in  Tkree- 
Dimeniional  Cylindrical  Geometry,"  Journal  of  Computational  Pkyiiet,  Vol.  78,  pp.  64-78,  1988.  Alio,  Strikwerda,  J.  C. 
and  Nagel,  Y  M.,  "A  Numemeal  Stud,  of  Flow  in  Spinning  and  Con, ny  Cylinder,,"  CRDC-SP-86007,  Proceeding,  of  the  1985 
Scientific  Conference  on  Chemical  Defente  Retearck,  Aberdeen  Proving  Ground,  Maryland,  April  1986. 

C--  'V“„'Ca'  ■'*,  J -  Fl*td  Dynamic,  Method,  for  Low  Reynold,  Number  Precctmg/ Spinning  Incompret.ible 

,ou#;  A  RBRL’MR*03657l  r  .S.  Arm y  Ballilite  ReaearzK  Laboratory,  Aberdeen  Proving  Grannd,  Maryland,  April  1968  (AD 
A  1 9 38 9 1  J, 
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1.  Cylindrical  Containers. 


Strikwerda  and  co-workers8  have  developed  a  steady-state,  finite-difference,  incom¬ 
pressible  Navier-Stokes  method  for  fully-filled  cylindrical  containers.  The  code  employs  an 
implicit,  iterative,  finite-difference  method  based  on  modified  line  successive-over- relaxation 
and  a  pressure  update  from  the  gradient  of  the  velocity  field.  This  code  (denoted  UWISC 
for  University  of  Wisconsin)  employs  a  non-uniform  grid  to  resolve  the  velocity  and  pres¬ 
sure  gradients  near  the  cylinder  walls,  central  finite-differences  in  the  radial  and  axial 
directions,  and  pseudo-spectral  differencing  to  represent  the  azimuthal  dependence.  Com¬ 
parison  of  UWISC  with  other  finite-difference  codes,  finite-element  codes,  and  analytical 
methods  is  reported  in  Reference  9. 

A  spatial  eigenvalue  method  has  been  developed  by  Hall,  Sedney,  and  Gerber10  (de¬ 
noted  HSG)  for  fully-filled  cylindrical  containers  and  constant- amplitude  coning  motion. 

Murphy  et.al.11  extended  this  method  to  partial-fill,  fully-filled  with  central  rod,  two 
immiscible  liquids,  and  slow  changes  to  the  amplitude  of  the  coning  motion.  The  Navier- 
Stokes  equations  are  written  in  an  inertial  reference  frame.  They  are  then  reduced  to  a 
set  of  linear  partial  differential  equations  based  on  small  angle  coning  motion,  and  linear 
departures  from  solid  body  rotation  are  considered.  No  slip  boundary  conditions  arc  im¬ 
posed  at  the  cylinder  walls  which  makes  all  flow  variables  proportional  to  where  / 

is  the  component  of  the  projectile  spin  along  the  symmetry  axis  and  t  is  time.  A  particular 
solution  was  employed  which  satisfies  the  axis  and  lateral  wall  boundary  conditions  but 
not  endwall  conditions.  The  eigenvalue  problem  is  obtained  using  separation  of  variables 
from  which  an  infinite  sequence  of  complex  eigenvalues  is  generated.  The  eigenvalues  are 
determined  by  an  iterative  process  for  which  sufficiently  accurate  initial  estimates  are  re¬ 
quired  for  convergence.  The  flow  variables  are  expressed  as  eigenfunction  expansions  with 
the  coefficients  determined  by  satisfying  the  endwall  boundary  conditions.  A  least  squares 
and  collocation  method  have  been  used  for  this  purpose.  Comparisons  of  measured  liq¬ 
uid  moment  coefficients  with  the  spatial  eigenvalue  code  (HSG)  and  UWISC  results9  have 
shown  the  consistency  of  both  methods.  However,  since  spatial  eigenvalue  methods  yield 
results  in  significantly  less  computer  run  time,  they  are  perhaps  the  preferred  scheme  for 
cylindrical  containers. 


2.  Non-cylindrical  Containers  -  High  Reynolds  Number  Flow. 

A  theory  for  non-cylindrical,  axisymmetric,  liquid-filled  containers  has  been  presented 
by  Wedemeyer.12  This  theory  is  derived  from  the  linearized,  invjscid  Navier-Stokes 
equations,  thus  Re  =  oo.  Using  this  theory,  the  approximate  e>gen-frequencies  of  the 
liquid  oscillations  can  be  computed.  Liquid  moments  occur  whenever  any  of  the  eigen- 
frequencies  fall  within  a  certain  bandwidth  of  the  fast  yaw  (precessional)  frequency  of  the 
projectile.  Liquid-filled,  non-cylindrical  cavities  can  be  examined  where  the  radius  of  the 

I0//a //,  P.,  Sedney,  R.,  and  Gerber,  N.,  “Fluid  Motion  in  Spinning,  Contnt  Cylinder  via  Spatial  Eigenfunction  Expansions, n 
A  RBRL-  TR'tSI 3,  US.  Army  Ballistic  Research  Laboratory,  Aberdeen  Proving  Ground,  Maryland,  August  1967. 

11  Murphy,  C.H.,  Bradley,  J.W.,  and  Mermagen,  W.H.,  "Liquid  Side  Moment  Exerted  by  a  Spinning,  Coning,  Highly  Viscous 
Liquid  Payload ,n  U.S.  Army  Ballistic  Research  Laboratory,  Aberdeen  proving  Ground,  Maryland,  report  in  preparation. 

12  Wedemeyer,  E.H.,  u Dynamics  of  Liquid  filled  Shell :  Non-Cylindncal  Cavity/'  BRL  Report  iVo.  1326,  U.S.  Army  Ballistic 

Research  Laboratory,  Aberdeen  proving  Ground,  Maryland,  1966. 
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cavity  is  a  slowly  varying  function  of  the  distance  along  the  cavity  axis,  |da/di|  1. 
For  cylindrical  cavities,  the  inviscid  eigen-frequencies  are  computed  exactly  and  given  in 
Stewaxtson  tables.13  For  small  deviations  from  cylindrical  shape,  small  changes  of  the 
eigen-frequencies  must  be  expected.  Karpov14  showed  experimentally  that  rounding  the 
corners  of  the  cylindrical  cavity  produced  very  little  effect  on  the  range  of  instability. 
However,  considerable  change  resulted  from  modifications  such  as  the  conical  reduction  of 
one  or  both  ends  of  the  cavity.  Wedemeyer’s  theory  agreed  favorably  with  Karpov’s  data. 


3.  Non-cylindrical  Containers  -  Low  Reynolds  Number  Flow. 

Written  at  the  U  S.  Army  Ballistic  Research  Laboratory  (BRL),  the  present  CFD 
code  is  a  modified  version  of  UWISC  (denoted  JWISC/BRL).  Like  the  UWISC  code,  the 
UWISC/BRL  code  uses  a  finite-difference  method  for  the  non-linear,  fully  viscous  Navier- 
Stokes  equations.  The  computational  algorithm  is  not  efficient  for  tlow  Reynolds  numbers 
above  300  (see  Section  IV. 2).  The  contained  flowfield  and  the  liquid-induced  moments  are 
computed  for  axisymmetric  containers  of  arbitrary  shape.  One  application  is  for  cylinders 
with  rounded  endcaps  where  |da/dr|  >  0.  The  endcap  is  a  circular  arc  where  the  cap 
height  above  the  flat  endwall,  at  the  axis  of  the  cylinder,  is  e/a.  For  a  hemispherical 
endcap,  e/a  =  1  (see  Figure  1). 


4.  Liquid  Roll  and  Side  Moment  Coefficients. 

Dimensional  analysis  and  linear  theories  indicate  that  the  liquid  moment  coeffi¬ 
cients  (roll  and  yaw  or  side)  will  depend  upon  the  following  dimensionless  parameters: 

Linear  Liquid  Moment  Coeff.  =  F[Re,  c/a,  r,  Kc/4>\  (1) 

Given  the  use  of  a  steady-state  CFD  code  that  retains  the  nonlinear  terms,  then 
the  liquid  moment  coefficient  will  also  depend  upon  the  precession  angle  ( Kc  =  sinac). 
However,  the  UWISC  steady-state  code  requires  Kc  =  0.  Hence,  this  code  would  vield  a 
dependence  as  follows: 

Steady  —  State  Nonlinear  Liquid  Moment  Coeff.  =  F(Re,  c/a,  r,  A'c]  (2) 

The  case  for  low  Re  should  also  follow  this  formulation  and  will  be  examined  using  the 
dimensionless  parameters  as  guides.  It  is  highly  possible  that  the  liquid  moment  coefficients 
are  linearly  related  to  ac  for  ac  <  20  degrees9.  If  this  is  the  case,  then  two  of  the  remaining 
three  parameters  can  be  held  constant,  while  the  behavior  of  the  liquid  moment  coefficient 
upon  the  third  parameter  can  be  explicitly  shown. 


•3  Suwariion.  A'.,  "On  The  Stability  0/  a  Sp tnninf  Top  Containing  Lignid,"  Journal  oj  Fluid  Mechanic),  Vol.  5,  Pi.  4. 
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Karpov  B.G.,  “Dynamic#  o)  Liqvid  Pilled  Shell.  Reaonance  in  Modified  Cylindrical  Catlli?#, "  BRL  Report  \o  13S2 
O  S.  Army  Balhttxe  Hue arc*  Laboratory,  Aberdeen  Proving  Ground,  Maryland.  Avgntt  <966. 
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Murphy15  suggested  the  use  of  roll  and  side  moment  coefficients  for  small,  fixed 
precession  angles,  ac,  as  defined  below: 

Roll  Moment  =  m/a2^>2[CLRMo  +  rAT2CLRM]  (3) 

Transverse  Moment  =  m/a24>2r[CLSM  +  * C\j,u\Kce'li‘c  (4) 


where, 

mi 

a 

4> 

T 

Ci.RM 

Clrm0 

Clsm 

Clim 

Ac 

<Pc 


is  the  mass  of  liquid  in  a  fully-filled  cylindrical  container  (2 irpia7c) 
is  the  maximum  radius  of  the  container 
is  the  spin  rate  of  the  container  in  the  inertial  frame 
is  the  ratio  of  coning  rate  to  spin,  <t>c/4> 

is  the  steady-state  liquid  roll  moment  coefficient  due  to  coning  motion 

is  the  liquid  roil  moment  coefficient  due  to  transient  liquid  effects  (i.e.  spinup) 

is  the  liquid  side  moment  coefficient 

is  the  liquid  in-plane  moment  coefficient 

is  sinac,  where  ac  is  the  precession  angle 

is  the  phase  angle  of  the  coiling  motion 


Further,  Reference  15  gives  a  relationship  between  the  moment  coefficients  for  the  lin¬ 
earized,  viscous  Navier-Stokes  equations.  Hence,  for  small  precession  angles  and  indepen¬ 
dent  of  the  Reynolds  number, 

Clrm  =  -Clsm  (5) 

For  arbitrary  precession  angles,  Rosenblat16  has  shown  that, 

Clrm  =  -CLSMtanac  (6) 


III.  COMPUTATIONAL  APPROACH 


In  tliis  section  the  implicit  finite-difference  method  is  described.  The  approach  fol¬ 
lows  that  of  Strikwerda8  for  a  cylindrical  container.  In  certain  aspects  of  the  method 
basic  changes  were  required  for  a  non-cylindrical  container  (e.g.  computational  grid,  grid 
stretching,  boundary  conditions).  In  particular  the  numerical  code  was  originally  written 
assuming  cylindrical  geometry  throughout.  Certain  other  aspects  of  the  method  are  ap¬ 
plicable  to  a  nou-cylindrical,  axisymmetric  container  (e.g.  psuedo-spectral  differencing). 
These  differences  and  similarities  in  the  method  are  highlighted  in  the  following  subsec¬ 
tions. 

1,1  Murphy,  C.  H.,  “A  Relation  Between  Liquid  Roll  Moment  and  Liquid  Side  Moment,’’  AIAA  Journal  of  Guidance,  Control 
and  Dynarr.ici,  VoL  8,  No.  t,  pp.  i87-t88,  March- April  1981.  Alto  ARBRL- MR-CSSdl,  U.5.  Arm y  tjaththc  Retearch 
Laboratory,  Aberdeen  Pronin g  'round,  Maryland,  April  I9id-  (AD  AldOSSd) 

,e  Roaenblat,  5.,  Gooding,  A.,  and  Englema n,  M.  S,,  “Finite  Element  Calculation i  of  Viaeo'ltetic  Flu  ’  F‘c  v  in  a  Spinning 
and  Nutating  Cylinder,”  CRDEC-CR-47011 ,  U.S.  Army  Chemical  Retearch,  Development  and  Engineering  Cenitr,  Aberdeen 
Proving  Ground,  Maryland,  December  lS8t. 
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1.  Equations  of  Motion 


Consider  a  fluid-filled  container  which  is  spinning  about  its  geometric  axis  with  spin- 
rate  ft,  and  about  the  coning  axis  with  spinrate  4>c  (see  Figure  1).  The  coning  axis  is  fixed 
in  the  inertial  (Earth-fixed)  reference  frame.  The  container  axis  is  inclined  to  the  coning 
axis  by  a  fixed  angle,  ac.  The  fluid  velocity  V  and  pressure  p  in  the  container  are  governed 
by  the  incompressible  Navier-Stokes  equations  for  constant  viscosity.  These  equations  in 
an  inertial,  cartesian  coordinate  system  are  given  by, 

^=-Vp+MV5?  (7) 


V-  V  =  0 


(8) 


These  equations  can  be  transformed  to  a  non-inertial  reference  frame  so  that  the  fluid 
motion  becomes  steady  and  thus  time  derivatives  are  zero,  d/dt  =  0.  This  reference  frame 
rotates  with  angular  velocity  <pc  about  the  coning  axis  and  was  also  used  by  Strikwerda.8 
Reference  8  gives  the  details  of  the  transformation.  In  addition,  the  cylindrical  coordinate 
system  is  used  in  which  ihe  radial,  azimuthal,  and  axial  directions  are  denoted  r,  <f>,  z  with 
velocity  components  it,  v,  w.  The  cylindrical  coordinate  system  is  also  used  in  the  endcaps 
where  the  grid  lines  are  not  aligned  with  the  r,  z  directions. 


The  equations  and  variables  are  non-dimensionalized  by  the  maximum  radius  of  the 
container  (a)  and  with  the  inertial  spinrate  (^  =  ft,  +  <j>ccosac)  and  the  container  radius, 
as  the  velocity  scale  (a<p).  Thus,  the  Reynolds  number  and  the  precessional  frequency  are 
defined  as, 


(9) 


r  = 


k 

<p 


(10) 


where  the  kinematic  viscosity  of  the  fluid  is  v  =  p!  pi- 


The  solid-body  rotation  is  subtracted  from  the  velocity,  and  the  pressure  is  redefined 
such  that, 

V  —  Krompuled  "b  I’to  lid-body  "b  T  X  P 


r2  r2 

P  =  Peomputed  +  y  4-  r3rcosac  +  —  [(rcos</>cosac  +  zsinac)7  -I-  rasin2ac] 
The  resultant  equations  are  given  by  Strikwerda8  as, 

1  dp 

A(ll)  ticonv  Ucorl  r ,  —  0 

Re  Or 

1  .  ,  ,  1  dp  n 

— A(u)  -  tw  -  vcor,  -  -  0 

^  Q 

—  A(w)  -  u>conv  -  u>cat\  ~  fa  ~  -2rsinacrcos</> 


(11) 

(12) 

(13) 

(14) 

(15) 
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where, 


1 3(ru)  1  dv  dw 
r  dr  '>r  r  d4>  dz 


du  v  du 

du  du 

Uconv  — 

udr  r  d<t>+Wdz  d4> 

r 

dv  v  dv 

dv  dv 

uv 

Vconv  — 

Ufr  +  ~rdi  +  Wd~z  +  d6- 

r 

dw  v  dw 

dw  dw 

tUeoav  = 

u—  +  “-5—  +u>—  4-  -r- 
dr  r  d<fr  dz  d<p 

«corf  =  2u>T8inac8in<^  —  2o(tcosotc  4-  1) 
Ujori  =  2u>Tsinaccos<£  +  2u(rcoso'c  4- 1) 
u>c«rl  —  -2rsinac(usin^  +  vcos<^) 


(16) 


1  d_ 

(  du\ 

1  dPu 

d*u 

2 

dv 

u 

A(u)  = 

r  dr 

V*) 

+  r3  d<f> 3  + 

dz 3 

r3 

d4> 

r3 

13 

(  dv\ 

1  d*v 

d3v 

2 

du 

V 

A(u)  = 

r  dr 

\dr) 

+  r3  d<P  + 

d?  + 

r3 

d</> 

'  r* 

A(m)  = 

1  d_ 

(  9w\ 

1  d*w 

d3w 

r  dr 

+  r3  d<t> 3  4 

'!& 

2.  Computational  Grid 


The  equations  of  motion  are  solved  on  a  discrete  grid  that  covers  the  entire  flow 
domain.  Figure  2  shows  a  typical  grid  for  a  cylinder  of  aspect  ratio  c/a  =  1.486  with 
rounded  endcaps  of  height,  e/a  =  1.  Only  the  portion  of  the  container  bounded  by  0  < 
r/a  <  l,<j>  =  0°,0  <  z/a  <  (c/a  +  e/a)  is  displayed  in  Figure  2.  In  this  example  there  are 
31  grid  points  in  the  r-direction,  81  grid  points  in  the  z-direction  and  19  grid  points  in 
the  endcaps  (i.e.  5°  spacing).  The  particular  grid  configuration  was  chosen  for  the  endcap 
region  to  insure  that  grid  cells  are  four-sided  (in  the  4>  —  constant  plane)  except  at  the 
axis,  r/a  =  0.  The  axis  grid  line  is  a  boundary,  along  which  boundary  conditions  are 
prescribed,  and  is  not  part  of  the  computational  domain.  Thus  multiple  grid  points  at  this 
location  are  inconsequential.  Grid  stretching  is  used  to  insure  adequate  grid  spacing  at 
the  walls  so  that  boundary  layers  and  viscous  shear  stresses  are  accurately  resolved.  This 
is  accomplished  using  tb*>  following  transformations,17 


+!)/(<* -1))>-1} 

1  +  [(<*+!)/(<*- 1))» 


(17) 


!’.V«»ca,  M.J.,  “Namarical  Simulation  oj  Unataady  Incomymtible  Flow  i»  t  Parlitllf.Filltd  Rolaimy  Cylindtr,”  ARBRL- 
TR-0t9ll,  U.S.  Arm)  Balliatie  Rtseartk  Laboratory,  Aktrdetn  Protiny  Ground,  Maryland,  Jam  1986  (AD  AI9600I). 
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(18) 


.  (e/a)(l-t  +  (l  +  t)|(t+  !)/(>-  Dr-1) 

Z/a  "  1  +  ((6  +  l)/(6  -  l))--1 

where  A  and  r)  are  uniformly  spaced  coordinates  in  the  r  and  z  directions,  respecively.  The 
parameters  d  <  1  and  b  <  1  axe  chosen  to  achieve  the  desired  clustering,  with  very  fine 
boundary  spacing  for  d  1  and  b  <£  1. 


3.  Computational  Algorithm 


The  numerical  method  used  to  solve  the  equations  of  motion  is  based  on  the  finite 
difference  scheme  of  Strikwerda18  and  a  psuedo-spectral  method.  Finite  differences  are 
used  to  approximate  derivatives  in  the  radial  (r)  and  axial  (r)  directions  and  the  psuedo- 
spectral  method  is  used  to  approximate  derivatives  in  the  azimuthal  direction  (<£).  Use  of 
the  psuedo-spectral  method  is  possible  since  the  flow  variables  vary  in  a  periodic  fashion  in 
the  4>  direction.  Since  the  non- cylindrical  geometry  under  consideration  in  this  work  is  also 
axisymmetric,  the  periodicity  in  4>  is  preserved  in  the  endcaps.  This  results  in  substantial 
savings  in  computer  memory  and  execution  time  requirements. 


Strikwerda18  has  shown  that  the  use  of  “regularized”  finite  differences,  instead  of 
central  difference  formula,  allows  the  equations  to  be  solved  on  a  highly  stretched  grid 
without  loss  in  accuracy.  As  examples  of  the  regularized  differences,  dp/dz  and  dw/dr  are 
approximated  by, 


dp 

_  Pk+l  ~  Pk~  1 

Pk+ 2  -  3p*+i  +  3p*  -  pk-\ 

dz 

2Ar 

6Az 

dw 

U>.+1  -  U><-.| 

ttfj+i  -  3 Wi  +  3u>i_i  -  u),_a 

2Ar 

6Ar 

(19) 

(20) 

where  t  and  k  represent  the  radial  and  axial  indicies,  respectively,  with  t  =  1  to  I  and 
k  =  1  to  K.  Grid  spacings  are  represented  by  Ar  and  A z.  Using  the  psuedo-spectral 
method,  approximations  to  the  pressure  derivatives,  for  example,  with  respect  to  the  an¬ 
gular  variable,  <p,  are  written  as8, 


(Jfi)' 

-  YL  na>,n,kCosn<t>j  -  nbi'^ksinrufij 

n=0 


(21) 


where  j  represents  the  azimuthal  indicie  with  j  —  1  to  J .  The  coefficients  are  given  by, 


2  J~l 

=  7  Y2  PuA8inn<^j 
J  ;=0 

2  *^-1 

bi,n,k  =  ~Z  ^  Pij,kCOSn<t>j 

J  j=0 


The  boundary  conditions  on  the  container  walls  follow  from  the  no-slip  condition  for 
velocity  and  extrapolation  from  the  interior  flow  for  the  pressure. 

_ u /■■),*  =  u»,>.i  =  “o.K  =  0 

‘‘Striku/erJa,  J.C.,  “/’inile  Difference  Meikoie  Jot  ikt  Stoka  and  Nantr-Slokei  Eqaatiom,”  SIAM  Journal  oj  Scientific 
and  Statulical  Comparing,  Vol.  S,  No.  March  1964- 
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VlJ.k  ~  vij,l  =  vij'K  =  0 

=  U =  Wij,K  =  0 

Pl.i.k  =  3PI-Ij,k  ~  3p/-2j,Jk  +P/-3J,* 
P« j, i  =  3p,j,a  -  3pi,j,3  +  pij.4 
P.j.tf  =  ZpiJ,K-l  -  Zpi,:,K-2  +  Pij,K-Z 


The  boundary  conditions  on  the  container  axis  (z-axis)  are  obtained  using  interpolation 
from  neighboring  grid  points,  for  example, 

1  1 4 


J.k  =  o  I  T 


1  /  4 


Plj,*  o 


Variables  at  the  grid  points  on  the  container  axis  at  the  interface  between  the  cylinder  and 
the  endcaps  (r  =  0,  s  —  |c/a|)  are  assigned  equal  values.  These  values  are  obtained  from 
interpolation  as  above. 

The  system  of  non-linear  equations  is  solved  using  a  modified  line  successive-over¬ 
relaxation  method  (LSOR).  The  method  is  described  by  Strikwerda19  for  the  case  of 
Linear  finite  difference  equations.  Because  the  psuedo-spectral  method  is  used  in  the  az¬ 
imuthal  direction,  line  relaxation  is  utilized.  In  this  case  each  line  contains  the  grid  points 
in  the  azimuthal  direction  for  each  value  of  the  radial  and  axial  coordinates.  This  approach 
is  unchanged  for  the  non-cylindrical  axisymmetric  geometry.  The  coupling  between  the 
velocity  components  and  the  natural  periodicity  of  the  azimuthal  coordinate  leads  to  a  pe¬ 
riodic  system  of  equations.  Thus  the  system  of  equations  that  must  be  solved  to  determine 
the  velocity  updates  for  each  radial  and  axial  grid  location  is  a  block  tridiagonal  periodic 
linear  system.  After  the  velocity  has  been  updated  by  one  pass  of  the  LSOR,  the  pressure 
is  updated  using, 

pm+ 1  =  pm  -  0  (V  •  Vm+1)  (22) 

where  m  is^tfce  iteration  index  and  /?  is  a  parameter.  Thus  V  •  V  =  0  is  not  solved  but 
rather  V  •  V  =  6  where  6  is  on  the  order  of  the  truncation  error  for  a  second-order  scheme. 


4.  Calculation  of  Liquid-Induced  Moments. 

Conservation  of  angular  momentum  for  the  steady  flow  in  a  control  volume  V  with 
surface  S  rotating  with  constant  angular  rate  about  a  fixed  axis  requires, 

M  /s(r-x/)dS 

=  Jv?x  (20  X  v)  pdV  +|vrx(fix(Qx  f)j  pdV  +  J  (r  x  9)  p9  •  dS  (23) 

'* Stnkmerd*,  J.C.,  “An  Iterntivt  Method  for  Solvmi  Finite  Difference  Affrotimtlione  to  the  Stokei  Efnttiont,"  SIAM 
Jonmol  of  Nnme ncal  Antlfm,  Voi  tl,  No.  3,  Jnne,  1934, 


8 


where  the  velocity  V  is  measured  relative  to  the  non-inertial  reference  frame.  In  this 
equation,  M  and  F  are  the  resultant  moment  on  the  control  volume  and  the  stress  acting 
on  the  walls  of  the  fully-filled  container,  respectively.  The  last  term  on  the  right-hand 
side  of  the  equation  vanishes  since  V  ■  dS  =  0  on  all  boundaries.  Thus  the  moment  can 
be  calculated  using  the  “surface  integral”  approach  or  the  “volume  integral”  approach. 
The  moment  can  be  expressed  in  terms  of  the  cartesian  components,  M  =  (M*,  Mv,  A/,), 
which  denote  the  yaw,  pitch  and  roll  moments,  respectively.  The  yaw  and  roll  components 
correspond  to  the  side  and  roll  moment  coefficients,  Clsm  and  Clrm-  Using  the  “volume 
integral”  approach,  Herbert20  and  Rosenblat16  have  shown  that  the  non-dimensional  yaw 
and  roll  moments  (non-dimensionalized  by  p/a5^2)  for  a  cylindrical  control  volume  can  be 
expressed  as, 


Mr  = 


2cosac 


J  J  J  (wr2 costfr'j  drdfidz 


M2  =  —Mxtsmac 


thus, 


Clsm  =  -  /  /  /  (u>r2cosd)  drd</>d* 

ir(c/ a)tanac  J-v  Jo  Jo  v  > 


(24) 

(25) 

(26) 


Clrm  =  -CLSMtanoc  (27) 

where  77  =  c/a  and  rmM  =  1.  Murphy11  has  also  derived  and  utilized  the  “volume  integral” 
in  the  spatial-eigenvalue  method  for  highly  viscous  liquids.  For  a  control  volume  consisting 
of  a  cylinder  with  rounded  endcaps,  Equations  26  and  27  are  used  with  rj  =  c/a  +  e/a  and 

t'xn&x  —  ^m*x(^)* 


IV.  RESULTS 

In  this  section  the  effects  of  rounded  endcaps,  versus  flat  endwalls,  on  the  liquid- 
induced  side  moment  coefficient,  Clsm5  is  examined  for  a  fully-filled  cylindrical  container. 
Results  are  shown  for  the  original  cylinder-only  version  of  the  finite-difference  code,  denoted 
UWISC,  and  the  Hall-Sedney-Gerber  spatial  eigenvalue  code,  denoted  HSG,  for  cylindrical 
geometries.  For  non-cylindncal  geometries,  results  from  the  generalized  axisymmetric 
geometry  version  of  the  finite-difference  code,  denoted  UWISC/BRL,  are  presented. 

Moments  were  computed  using  both  the  “surface  integral”  and  the  “volume  integral” 
approach  (see  Section  III. 4).  These  methods  yielded  equivalent  numerical  results  for  all 
cases  presented  in  this  report.  However,  the  “volume  integral”  approach  requires  less 
computational  effort. 

Because  the  HSG  code  is  formulated  assuming  a  small  coning  angle,  all  of  the  UWISC 
and  UWISC/BRL  computations  were  performed  for  ac  =  2°.  All  UWISC  computations 
presented  in  this  section  were  performed  using  the  same  grid  dimensions  with  variable  grid 
wall  spacing.  All  HSG  computations  were  performed  using  15  to  40  eigenvalues  for  Re  21.5 
to  300.  In  some  cases  a  cubic  spline  was  used  to  join  discrete  data  points  in  the  following 
figures. 

30 Herbert.  T.,  “On  the  Vi,co*t  Roll  Moment  in  a  Spinninp  ini  Nutting  Cylinder, "  CRDC-SP.S600T,  Proceed, n/>  oj  the 
1 9S4  Scientific  Conference  on  Chemical  Defense  Research,  Aberdeen  Proving  Ground,  Maryland,  April  1985, 
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1.  Variation  of  Clsm  with  c/a. 


Figure  3  shows  the  variation  of  Clsm  with  endcap  height,  e/a,  for  a  cylinder  of  aspect 
ratio,  c/a  =  1,486.  The  height  of  the  endcap  is  measured  at  the  cylinder  axis  (see  Figure 
1).  For  e/a  =  0  the  container  is  a  cylinder  with  flat  endwalls,  and  for  e/a  =  1  the  endcaps 
are  hemispher  cal.  The  same  value  of  Clsm  for  the  cylindrical  geometry  is  produced  by 
the  UWISC/BRL  and  HSG  codes.  The  value  of  Clsm  decreases  with  increasing  endcap 
height,  and  reaches  a  minimum  for  e/a  —  1.  This  suggests  that  the  liquid  side  moment  can 
be  reduced  by  as  much  as  20%  by  rounding  the  endwalls  of  a  cylinder  with  aspect  ratio 
c/a  =  1.486  at  Re  =  50.  The  value  of  the  coning  frequency  used  in  this  figure  (r  =  .0469)  is 
typical  of  ballistic  projectiles  with  liquid  payloads.  Since  an  endcap  height  of  1  represents 
the  limiting  case,  all  subsequent  calculations  are  for  e/a  =  1. 


2.  Variation  of  Clsm  with  Re. 

Figure  4  shows  the  variation  of  Clsm  with  Reynolds  number  for  endcap  heights  of 
0  and  1,  and  a  cylinder  of  aspect  ratio  of  1.486.  For  the  cylindrical  geometry,  both  the 
UWISC/BRL  and  HSG  codes  show  the  same  variation  of  liquid  side  moment  with  Re, 
including  a  maximum,  or  resonance,  at  about  Re  =  75.  Results  using  the  UWISC/BRL 
code  for  the  same  cylinder  with  hemispherical  endcaps  show  a  resonant  condition  at  about 
Re  =  50.  In  addition,  the  value  of  Clsm  at  resonance  is  about  22%  smaller  for  the  cylinder 
with  rounded  endcaps.  For  Re  =  21.5  the  liquid  side  moment  is  actually  12%  smaller  for 
the  cylinder  with  flat  endwalls.  However,  for  Re  greater  than  the  resonance  value,  the 
addition  of  rounded  endcaps  has  reduced  the  Clsm  by  as  much  as  55%. 

As  noted  in  a  previous  study9  the  UWISC  code  becomes  computationally  inefficient 
for  for  Re  >  300.  Figure  5  shows  the  variation  of  computer  (CRAY-XMP/48)  CPU  time 
with  Reynolds  number.  In  each  case  the  solution  for  a  particular  Re  was  used  as  a  starting 
solution  for  the  subsequent  Re  solution.  Since  the  computer  time  required  increases  in  an 
almost  linear  fashion  with  Re  when  scratch  starts  are  used,  the  restart  method  results  in  a 
substantial  savings  in  computer  run  time.  For  Re  <  100,  solutions  can  be  achieved  in  2  to 
4  hours  using  about  30,000  grid  points  for  the  cylindrical  geometry,  and  7000  additional 
grid  points  for  the  rounded  endcaps.  These  values  are  considered  more  than  adequate  for 
c/a  =  1.486.  In  the  previous  study9  with  cylindrical  containers,  computer  run  times  were 
slightly  less  than  2  hours  for  the  same  value  of  c/a  since  a  smaller  grid  was  used.  The 
Re=100  solution  for  the  cylindrical  container  required  slightly  more  computer  time  than 
the  cylinder  with  rounded  endcaps.  The  data  from  Figure  4  indicate  that  the  maximum 
Clsm  also  occurs  at  about  Re=100.  This  indicates  that  rounding  the  cylinder  endcaps  may 
promote  faster  solution  convergence.  For  larger  values  of  Re,  the  computer  run  time  rapidly 
increases  beyond  5  hours  and  approaches  12  hours  for  Re  =  300.  Some  reduction  in  run  time 
(<5%)  can  be  realized  by  changing  the  value  of  $  in  Equation  22  (i.e.  a  converged  solution 
would  require  fewer  iterations)  or  relaxing  the  convergence  tolerance.  Part ially- converged 
solutions  (i.e.  using  relaxed  convergence  tolerance)  for  Reynolds  Numbers  around  10  can 
require  <2  CPU  hours  using  a  mini-computer  (e.g.  VAX  8600)  and  be  used  for  preliminary 
design  applications.  This  technique  results  in  severe  degradation  of  the  accuracy  of  Clsm 
values  for  Re  greater  than  about  20. 
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3.  Variation  of  Clsm  with  r. 

Figure  6  shows  the  variation  of  Clsm  with  coning  frequency  for  Re  =  50,  endcap 
heights  of  0  and  1,  and  a  cylinder  of  aspect  ratio  of  1.486.  For  the  cylindrical  geometry, 
both  the  UWISC/BRL  and  HSG  codes  show  the  same  variation  of  liquid  side  moment 
with  r,  including  a  resonance  at  r  =  .5.  Results  using  the  UWISC/BRL  code  for  the  same 
cylinder  with  hemispherical  endcaps  show  a  resonant  condition  at  r  =  .65,  with  a  value  of 
Clsm  about  20%  larger  than  the  Clsm  for  a  cylinder.  For  r  greater  than  about  .35,  the 
liquid  side  moment  is  actually  50%  smaller  for  the  cylinder  with  flat  endwalls.  However, 
for  r  smaller  than  the  resonance  value,  the  addition  of  rounded  endcaps  has  reduced  the 
Clsm  by  as  much  as  19%. 

As  noted  in  a  previous  study9  the  UWISC  code  cannot  be  run  for  r  >  1,  whereas 
the  HSG  code  has  been  used  for  coning  frequencies  greater  than  the  spin  frequency.  This 
condition  can  be  resolved  by  reformulating  the  UWISC  code  and  dividing  the  governing 
equations  by  r. 


4.  Variation  of  Clsm  with  c/a. 

Figure  7  shows  the  variation  of  Clsm  with  cylinder  aspect  ratio  for  Re  =  50,  r  =  .0469, 
and  endcap  heights  of  0  and  1.  For  the  cylindrical  geometry,  both  the  UWISC/BRL  and 
HSG  codes  show  the  same  variation  of  liquid  side  moment  with  c/a,  including  a  resonance 
at  c/a  =  1.2.  Results  using  the  UWISC/BRL  code  for  the  same  cylinder  with  hemispherical 
endcaps  show  a  resonant  condition  at  c/a  =  .6,  at  about  the  same  value  of  Clsm.  For  c/a 
greater  than  1.0  the  liquid  side  moment  is  up  to  19%  smaller  for  the  cylinder  with  rounded 
endcaps.  For  c/a  greater  than  2  the  value  of  Clsm  is  only  2%  smaller  for  cylinders  with 
rounded  endcaps.  The  value  of  Clsm  approaches  0  as  c/a  approaches  0  for  both  the 
cylinder  with  flat  endwalls  and  the  cylinder  with  hemispherical  endcaps.  In  the  case  of  the 
hemispherical  endcaps,  for  c/a  =  0  the  container  is  a  sphere.  For  a  sphere,  the  integration 
of  pressures  over  the  symmetric  interior  geometry  will  not  produce  a  contribution  to  Clsm- 
Additionally,  viscous  shear  stresses  will  be  zero  since,  for  a  sphere,  steady  rotation  and 
precession  is  simply  rigid  body  rotation  about  a  different  axis. 


5.  Application  to  Flight  Tested  Payloads. 

Pope21  reported  analysis  of  flight  test  data  for  projectiles  with  highly  viscous  liquid 
payloads  (i.e.  low  Reynolds  Number)  contained  in  cylinders  with  large  c/a.  A  particular 
example  was  for  Re  =  45.2,  r  =  .123,  c/a  =  4.23,  and  ae  =  2°.  Table  1  lists  the  measured 
Clsm  33  well  as  computed  Clsm  using  the  HSG  and  UWISC/BRL  codes.  Both  HSG  and 
UWISC/BRL  show  a  Clsm  value  23%  smaller  than  measured  for  the  cylindrical  container. 
The  addition  of  rounded  endcaps  slightly  increases  the  moment  for  e/a  <  1  and  decreases 
the  moment  by  <1%  for  e/a  =  1.  The  small  change  in  Clsm  for  this  case  is  consistent 
with  the  results  of  Figure  7  for  large  aspect  ratio  containers.  In  this  case  “rounding” 

31  Pope,  R.L.,  “Further  Analyiit  of  Yaweonie  Data  From  Some  Liqaid  Payload  Projectilea,"  ARBRL-MR-033t9,  U.S,  Army 
Balliittc  Reeearth  laboratory,  Aberdeen  Proviny  Groand,  Maryland,  December  1933  (AD  At37tSS). 
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the  endcaps  does  not  substantially  reduec  the  liquid  moment  and  would  not  represent  an 
effective  solution  for  unstable  projectile  flights. 


Table  1.  Application  to  Flight  Tested  Payloads. 


t/a 

Datum 

Clsm 

0 

Flight  Data31 

.0550 

0 

HSG 

.0424 

0 

UWISC/BRL 

.0420 

.50 

UWISC/BRL 

.0424 

.75 

UWISC/BRL 

.0421 

1.0 

UWISC/BRL 

.0417 

V.  CONCLUSIONS 


The  three-dimensional,  steady,  laminar,  Navier-Stokes  equations  are  solved  using  an 
implicit  finite-difference  scheme  based  on  successive-over-relaxation.  The  method  of  Strik- 
werda,  for  fluid-filled  cylindrical  containers,  has  been  extended  to  non-cylindrical  axisym- 
metric  geometries.  The  liquid  side  moment  coefficient  for  a  cylinder  with  rounded  endcaps 
is  investigated  using  this  code  for  Reynolds  number,  10  <  Re  <  300.  Rounded  endcaps  can 
decrease  the  liquid-induced  moment  by  as  much  as  40%  for  some  values  of  the  Reynolds 
number,  coning  frequency,  and  cylinder  aspect  ratio.  . 

Further  modification  of  the  code  (e.g.  ellipsoidal  endcaps,  partial-fill,  etc.)  could  be 
accomplished.  However,  a  computational  fluid  dynamics  code  for  time-dependent,  three- 
dimensional,  incompressible  flow  in  arbitrary  geometries  ha>:  been  written  by  Chakravarthy 
et.al.23.  This  code  is  a  generalized  version  of  a  previous  code  ‘23'2A  that  was  used  for 
fully-filled  cylindrical  containers  in  steady  coning  motion  and  flow  Reynolds  Numbers  as 
high  as  500,000.  The  new  code  by  Chakravarthy  uses  the  same  numerical  scheme  as  the 
previous  code  but  includes  a  zonal  grid  scheme  to  faciliate  grid  generation  in  containers 
of  arbitrary  shape.  The  code  requires  modification  of  the  boundary  conditions,  reference 
frame,  and  equations  of  motion  in  order  to  compute  liquids  in  precessing/spinning  contain¬ 
ers.  However,  this  code  can  be  more  easily  modified  than  UWISC/BRL  for  free  surface 
boundary  conditions,  multiple  liquids,  time-dependent  container  motion  and  other  cases 
of  interest. 


23  Ptn,  D-i  Ckakravarthp,  S.R.,  “Unified  Formulation  /or  Ineomprtaaikle  Flowe,"  AIAA-89-0ltt,  Proceeding  of  ike 
t7tk  AIAA  Aeroepace  Science*  Meet  inf,  Reno,  NV,  Jtnvery  9. it,  1989. 

33  Kwak,  D.,  end  Chakravarthy,  S.R.,  “A  Tkret-Dimeneional  Incompreeeihle  Navier-Siokee  Flow  Solver  Uetny  Primitive 
Variablce, "  AIAA  Jonrnei  Vol.  ti,  No.  3,  Merck  1989,  pp.  390-398. 

34  Nueca,  M.J.,  end  Cooper,  G.R.,  “Compvtetionel  Fluid  Dpnemtce  Method /or  Vneteadp  IncomprttiiHe  Rotatiny  Flow  wtth 
/mptliite  Con  in/  Motion,’'  CRDEC-SP-8900S,  Proceeding  oj  the  1988  Scientific  Conference  on  Chemical  Dtjenae  Reieerch, 
Aherdetn  proviny  Ground,  Maryland,  1989. 
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Figure  1.  Cylinder  With  Rounded  Endcaps  -  Configuration  and  Nomenclature  (Hemi¬ 
spherical  Endcap  Shown) 
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Figure  2.  Typical  Computational  Grid  (Hemispherical  Endcap  Shown) 
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Re  =  50  t  =  .0469  c/a  =  1.486  ae  =  2deg 
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Figure  3.  Variation  of  Liquid-Induced  Side  Moment  Coefficient  with  Rounded  Endcap 
Height;  Re  =  50,  r  =  .0469,  c/a  =  1.486,  ac  =  2° 
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Figure  4.  Variation  of  Liquid-Induced  Side  Moment  Coefficient  with  Reynolds  Number; 
e/a  =  0  and  1,  r  =  .0469,  c/a  =  1.486,  ac  ~  2° 
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Figure  5.  Variation  of  CRAY-XMP/48  CPU  Time  with  Reynolds  Number;  e/a  =  0  and 
1,  r  =  .0469,  c/a  =  1.486,  ac  =  2° 
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Figure  6.  Variation  of  Liquid-Induced  Side  Moment  Coefficient  with  Coning  Frequency; 
c/a  =  0  and  1,  He  =  50,  c/a  =  1.486,  ac  =  2° 
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Figure  7.  Variation  of  Liquid-Induced  Side  Moment  Coefficient  with  Cylinder  Aspect 
Ratio;  e/a  -  0  and  1,  Re  =  50,  r  =  .0469,  ac  =  2° 


References 


1.  Murphy,  C.H.,  “Angular  Motion  of  a  Spinning  Projectile  with  a  Viscous  Liquid  Pay- 
load,”  ARBRL-MR-03194,  U.S.  Army  Ballistic  Research  Laboratory,  Aberdeen  Prov¬ 
ing  Ground,  Maryland,  August  1982  (AD  A118676).  Also,  AIAA  Journal  of  Guidance, 
Control,  and  Dynamics,  Vol.  6,  No.  4,  pp. 280-286,  July-August  1983. 

2.  Gerber,  N.  and  Sedney,  R.  “Moment  on  a  Liquid-Filled  Spinning  and  Nutating  Pro¬ 
jectile:  Solid  Body  Rotation,”  ARBRL-TR-02470,  U.S.  Army  Ballistic  Research  Lab¬ 
oratory,  Aberdeen  Proving  Ground,  Maryland,  February  1983  (AD  A125332). 

3.  Murphy,  C.  H.,  “Moment  Induced  by  a  Liquid  Payload  During  Spin-Up  Without  a 
Critical  Layer,”  ARBRL-TR-02581,  U.S.  Army  Ballistic  Research  Laboratory,  Ab¬ 
erdeen  Proving  Ground,  Maryland,  August  1984  (AD  A 145716).  Also,  AIAA  Journal 
of  Guidance,  Control,  and  Dynamics,  Vol.  8,  No.  3,  pp.354-359,  May-June  1985. 

4.  Gerber  N.,  “Liquid  Moment  on  a  Filled  Coning  Cylinder  During  Spin-Up:  Ad  Hoc 
Model,”  ARBRL-TR-2628,  U.S.  Army  Ballistic  Research  Laboratory,  Aberdeen  Prov¬ 
ing  Ground,  Maryland,  December  1984  (AD  150280). 

5.  D’Amico  W.  P.,  “Flight  Data  on  Liquid-Filled  Shell  for  Spin-Up  Instabilities,”  ARBRL- 
MR-03334,  U.S.  Army  Ballistic  Research  Laboratory,  Aberdeen  Proving  Ground, 
Maryland,  February  1984  (AD  139136).  Also,  AIAA  Paper  83-2143,  August  1983. 

6.  Nusca  M.  J.,  D’Amico  W.  P.,  and  Beims,  W.  G.,  “Pressure  Measurements  in  a  Rapidly 
Rotating  and  Coning,  Highly  Viscous  Fluid,"  ARBRL-MR-03325,  U.S.  Army  Ballistic 
Research  Laboratory,  Aberdeen  Proving  Ground,  Maryland,  November  1983  (AD 
A136824). 

7.  Hepner,  D.J.,  Kendall,  T.M.,  Davis,  B.S.,  and  Tenly,  W.Y.,  “Pressure  Measure¬ 
ments  in  a  Liquid-Filled  Cylinder  Using  A  Three- Degree-Of  Freedom  Flight  Simu¬ 
lator,”  ARBRL-MR-3560,  U.S.  Army  Ballistic  Research  Laboratory,  Aberdeen  Prov¬ 
ing  Ground,  Maryland,  December  1986.  Also,  AIAA  Paper  No.  86-2026,  1986  (AD 
A177872). 

8.  Strikwerda,  J.  C.,  and  Nagel,  Y.  M.,  “A  Numerical  Method  for  the  Incompressible 
Navier-Stokes  Equations  in  Three-Dimensional  Cylindrical  Geometry,” 

Journal  of  Computational  Physics,  Vol.  78,  pp.  64-78,  1988.  Also,  Strikwerda,  J.  C., 
and  Nagel,  Y.  M.,  “A  Numerical  Study  of  Flow  in  Spinning  and  Coning  Cylinders,” 
CRDC-SP-86007,  Proceedings  of  the  1985  Scientific  Conference  on  Chemical  Defense 
Research,  Aberdeen  Proving  Ground,  Maryland,  April  1986. 

9.  Nusca,  M.  J.,  “Computational  Fluid  Dynamics  Methods  for  Low  Reynolds  Num¬ 
ber  Precessing/Spinning  Incompressible  Flows,”  ARBRL-MR-03657,  U.S.  Army  Bal¬ 
listic  Research  Laboratory,  Aberdeen  Proving  Ground  Maryland,  April  1988  (AD 
A193891). 

10.  Hall,  P.,  Sedney,  R.,  and  Gerber,  N.,  “Fluid  Motion  in  Spinning,  Coning  Cylinder  via 
Spatial  Eigenfunction  Expansions,”  ARBRL-TR-2813,  U.S.  Army  Ballistic  Research 
Laboratory,  Aberdeen  Proving  Ground,  Maryland,  August  1987. 


20 


11.  Murphy,  C.H.,  Bradley,  J.W.,  and  Mermagen,  W.H.,  “Liquid  Side  Moment  Exerted 
by  a  Spinning,  Coning,  Highly  Viscous  Liquid  Payload,”  U.S.  Army  Ballistic  Research 
Laboratory,  Aberdeen  Proving  Ground,  Maryland,  report  in  preparation. 

12.  Wedemeyer,  E.H.,  “Dynamics  of  Liquid  Filled  Shell:  Non- Cylindrical  Cavity,”  BRL 
Report  No.  1326,  U.S.  Army  Ballistic  Research  Laboratory,  Aberdeen  Proving  Ground, 
Maryland,  August  1966. 

13.  Stewartson,  K.,  “On  The  Stability  of  a  Spinning  Top  Containing  Liquid,” 

Journal  of  Fluid  Mechanics.  Vol.  5,  Pt.  4,  1959. 

14.  Karpov,  B.G.,  “Dynamics  of  Liquid  Filled  Shell:  Resonance  in  Modified  Cylindrical 
Cavities,”  BRL  Report  No.  1332,  U.S.  Army  Ballistic  Research  Laboratory,  Aberdeen 
Proving  Ground,  Maryland,  August  1966. 

15.  Murphy,  C.  H.,  “A  Relation  Between  Liquid  Roll  Moment  and  Liquid  Side  Moment,” 
AIAA  Journal  of  Guidance,  Control  and  Dynamics,  Vol.  8,  No.  2,  pp.  287-288, 
March-April  1985.  Also,  ARBRL-MR-03347  U.S.  Army  Ballistic  Research  Labora¬ 
tory,  Aberdeen  Proving  Ground,  Maryland,  April  1984  (AD  A140658). 

16.  Rosenblat,  S.,  Gooding,  A.,  and  Engleman,  M.  S.,  “Finite  Element  Calculations  of 
Viscoelastic  Fluid  Flow  in  a  Spinning  and  Nutating  Cylinder,”  CRDEC-CR-87021, 
U.S.  Army  Chemical  Research,  Development  and  Engineering  Center,  Aberdeen  Prov¬ 
ing  Ground,  Maryland,  December  1986. 

17.  Nusca,  M.J.,  “Numerical  Simulation  of  Unsteady  Incompressible  Flow  in  a  Partially- 
Filled  Rotating  Cylinder,”  ARBRL-TR-02915,  U.S.  Army  Ballistic  Research  Labora¬ 
tory,  Aberdeen  Proving  Ground,  Maryland,  June  1988  (AD  A196001). 

18.  Strikwerda,  J.C.,  “Finite  Difference  Methods  for  the  Stokes  and  Navier-Stokes  Equa¬ 
tions,”  SIAM  Journal  of  Scientific  and  Statistical  Computing,  Vol.  5,  No.  1,  March 
1984. 

19.  Strikwerda,  J.C.,  “An  Iterative  Method  for  Solving  Finite  Difference  Approximations 
to  the  Stokes  Equations,”  SIAM  Journal  of  Numerical  Analysis,  Vol.  21,  No.  3,  June, 
1984. 

20.  Herbert,  T.,  “On  the  Viscous  Roll  Moment  in  a  Spinning  and  Nutating  Cylinder,” 
CRDC-SP-86007,  Proceedings  of  the  1984  Scientific  Conference  on  Chemical  Defense 
Research,  Aberdeen  Proving  Ground,  Maryland,  April  1985. 

21.  Pope,  R.L.,  “Further  Analysis  of  Yawsonde  Data  From  Some  Liquid  Payload  Pro¬ 
jectiles,”  ARBRL-MR-03329,  U.S.  Army  Ballistic  Research  Laboratory,  Aberdeen 
Proving  Ground,  Maryland,  December  1983  (AD  A137256). 

22.  Pan,  D.,  and  Chakravarthy,  S.R.,  “Unified  Formulation  for  Incompressible  Flows,” 
AIAA-89-0122,  Proceedings  of  the  27th  AIAA  Aerospace  Sciences  Meeting,  Reno, 
NV,  January  9-12,  1989. 

23.  Kwak,  D.,  and  Chakravarthy,  S.R.,  “A  Three-Dimensional  Incompressible  Navier- 
Stokes  Flow  Solver  Using  Primitive  Variables,”  AIAA  Journal.  Vol.  24,  No.  3,  March 
1986.  pp.  390-396. 


21 


24.  Nusca,  M.J.,  and  Cooper,  G.R.,  “Computational  Fluid  Dynamics  Method  for  Un¬ 
steady  Incompressible  Rotating  Flow  with  Impulsive  Coning  Motion,”  CRDC-SP- 
89005,  Proceedings  of  the  1988  Scientific  Conference  on  Chemical  Defense  Research, 
Aberdeen  Proving  Ground,  Maryland,  1989. 


22 


List  of  Symbols 


a  cylinder  maximum  radius 

c  cylinder-section  half-height 

c/a  cylinder-section  aspect  ratio 

Clrm  steady-state  liquid  roll  moment  coefficient  due  to  coning  motion 

Clrmo  liquid  roll  moment  coefficient  due  to  transient  liquid  spinup 

Clsm  liquid  side  moment  coefficient 

Clim  liquid  in-plane  moment  coefficient 

e/a  ratio  of  endcap  height  to  maximum  cylinder  radius 

F  shear  stress  on  the  container  surface 

i,j,k  grid  indicies  for  the  radial,  azimuthal,  and  axial  directions 

I,  J,  K  total  number  of  grid  points  in  the  radial,  azimuthal,  and  axial  directions 

Ke  sinac 

Kc  yaw  growth  rate 

M  angular  momentum 

mi  mass  of  liquid  in  a  fully-filled  cylindrical  container  (27rp/a2e) 
p  static  pressure 

r  radius  or  radial  direction 

r  position  vector 

Re  Reynolds  number,  a?</>/v 

S  control  surface 

t  time 

u,  v,  w  velocity  components  in  the  radial,  azimuthal,  and  axial  directions 

V  velocity  vector 

V  control  volume 

z  axial  direction 

Greek  Symbols 


<*c 

& 

M 

i/ 

Pi 

T 

<? 

4>c 

j> 

0c 

Q, 


coning  or  precession  angle 
parameter,  see  Equation  22 
molecular  viscosity  of  the  liquid 
kinematic  viscosity  of  the  liquid,  p/p 
density  of  the  liquid 
ratio  of  coning  rate  to  spin, 
azimuthal  angle  or  azimuthal  direction 
phtve  »ngle  of  the  coning  motion 

spin  rate  of  the  container  in  the  inertial  frame,  +  <^ccosac 
spin  rate  of  the  container  about  the  coning  axis 
spin  rate  of  the  container  about  the  geometric  axis 
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Intentionally  left  blank. 
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SMCCR-MU 
SMCCR-MSI 

Dir.  VLAMO 

ATTN:  AMSLC-VL-D 
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Space  Administration 
Ames  Research  Center 
ATTN:  Dr.  J.  Steger 
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1  Aerospace  Corporation 

Aero- Engineering  Subdivision 
ATTN:  Walter  F.  Reddall 
El  Segundo,  CA  90245 

1  Calspan  Corporation 
ATTN:  W.  Rae 
P.O.  Box  400 
Buffalo,  NY  14225 

1  Hughes  Aircraft 

ATTN:  Dr.  John  McIntyre 

Mail  Code  S41/B323 
P.O.  Box  92919 
Los  Angeles,  CA  90009 

1  Rockwell  International 
Science  Center 
ATTN:  Dr.  S.  Giakravarthy 
1049  Camino  Dos  Rios 
P.O.  Box  1085 
Thousand  Oaks,  CA  91360 

1  Arizona  State  University 
Department  of  Mechanical  and 
Energy  Systems  Engineering 
ATTN:  G.P.  Neitzel 
Tempe,  AZ  85281 

1  Illinois  Institute  of  Technology 
ATTN:  Mr.  Simon  Rosenblat 
3300  South  Federal 
Chicago,  IL  60616 

1  Massachusetts  Institute  of 
Technology 
ATTN:  H.  Greenspan 
77  Massachusetts  Avenue 
Cambridge,  MA  02139 
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1  North  Carolina  State  University 
Mechanical  and  Aerospace 
Engineering  Department 
ATTN:  F.F.  DeJamette 
Raleigh,  NC  27607 

1  Northwestern  University 
Dept,  of  Engineering  Science 
and  Applied  Mathematics 
ATTN:  Dr.  S.H.  Davia 
Evanston,  IL  60201 

1  Chio  State  University 

Dept,  of  Mechanical  Engineering 
ATTN:  Dr.  T.  Herbert 
Columbus,  OH  43221 

l  Renssalaer  Polytechnic  Institute 
Department  of  Math  Sciences 
Troy,  NY  12181 

1  University  of  Oolorado 

Department  of  Astro-Geophysics 
ATTN:  E.R.  Benton 
Boulder,  CO  80302 

1  University  of  Maryland 
ATTN:  J.D.  Aider  son 
College  Park,  MD  20740 

1  University  of  Maryland 

Baltimore  County 
Department  of  Mathematics 
ATTN:  Dr.  Y.M.  Lynn 
5401  Wilkens  Avenue 
Baltimore,  MD  21228 

2  Uhiv.  of  Southern  California 
Dept,  of  Aerospace  Engineering 
ATTN:  T.  Maxworthy 

P.  Weidman 

Los  Angeles,  CA  90007 

1  University  of  Wisconsin-Mad iscn 
Center  for  Mathematical  Sciences 
ATTN:  Dr.  John  Strikwerda 
610  Walnut  Street 
Madison,  WI  53706 

1  VPI  and  State  University 

Dept,  of  Aerospace  Engineering 
ATTN:  Tfech  Library 
Blacksburg,  VA  24061 


26 


AUTHOR'S  DISTRIBUTION 


No.  of 

Copies  Organization 

1  Mr.  Harold  Vaughn 
7709  Gladden  N.E. 
Albuquerque,  M4  87110 


Aberdeen  Proving  Ground 


Director,  USAMSAA 
ATTN:  AMXSY-RA,  R.  Scungio 

Conmander,  USATECOM 
ATTN:  AMSTE-TE-F,  W.  Vcmocil 
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W.  Dee 

C.  Hughes 
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